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We present measured dielectron production cross sections for Ca+Ca, C+C, 
He+Ca, and d+Ca reactions at 1.0 A-GeV. Statistical uncertainties and sys- 
tematic effects are smaller than in previous DLS nucleus-nucleus data. For 
pair mass M < 0.35 GeV/c 2 : 1) the Ca+Ca cross section is larger than the 
previous DLS measurement and current model results, 2) the mass spectra 
suggest large contributions from 7r° and r\ Dalitz decays, and 3) da/dM oc 
A P -A T . For M > 0.5 GeV/c 2 the Ca+Ca to C+C cross section ratio is 
significantly larger than the ratio of Ap-Ap values. 
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Dielectrons produced in heavy-ion collisions are attractive probes for studying dynam- 
ical properties of nucleus-nucleus interactions. The e + e~ pairs do not undergo significant 
rescattering in the reaction, thus the kinematics of the pairs retains information about their 
production. This is of particular interest if the e + e~ pairs are produced by processes, such 
as pion annihilation, that must occur during in the hot, dense phase of the collisions. Use 
of this probe has produced interesting results at both Bevalac Q and SPS energies. We 
present in this letter the latest measurements of dielectron production from the Dilepton 
Spectrometer (DLS) Collaboration in nucleus-nucleus reactions at a beam kinetic energy of 
1.0 A-GeV. 

The DLS collaboration has previously reported on dielectron production in several col- 
liding systems [l],3,f|. The first generation DLS data from p+Be, Ca+Ca, and Nb+Nb |]|| 
reactions provided the first observations of dielectrons produced at Bevalac energies. Early 
calculations suggested that such data could be dominated by contributions from 7r + 7r~ anni- 
hilation . Subsequent models of AA collisions in this energy regime |7|-[| calculated that 
e + e~ pairs of invariant mass below about 0.4 GeV/c 2 are produced primarily from conven- 
tional hadronic sources, such as pn bremsstrahlung and Dalitz decay processes (tt°, A, and 
rj), but that contributions from 7r + 7r~ annihilation were needed to explain the Ca+Ca data 
at higher pair masses. Models that focus on density induced changes in the p-meson mass 
provide alternative descriptions of the pair yield at the higher masses [HIE!- Within the 



limited statistics of the first generation DLS data, it was not possible to distinguish among 
the models that provided results for specific DLS measurements. 



After improvements to the DLS apparatus [13-15], a second generation of measurements 



was obtained: first from p+p and p+d reactions at a number of energies ||, and then from 
the Ca+Ca, C+C, He+Ca, and d+Ca reactions presented in this letter. Each of these data 
sets contain significantly more pairs than earlier DLS data. To increase our sensitivity to 
the effects of multiple hadronic interactions (e.g. 7r + 7r~ annihilation and multi-step reso- 
nance excitation), the nucleus-nucleus reactions were chosen to have different numbers of 
participant nucleons, but identical isospin and similar internal nuclear motion. 



A description of the DLS apparatus has been published fl3fl , and the analysis procedures 



used here are discussed in Ref. [jM]. The new procedures lead to a mass-independent mass 
resolution for these data of 10%, and smaller bin widths for M < 0.2 GeV/c 2 . No explicit 
requirement on impact parameter (e.g. multiplicity) is imposed on these data. The cross 
sections presented here are evaluated in a three-dimensional, binned array of invariant mass 
(M), transverse momentum (pt), and laboratory rapidity (y). The data are available from 
the authors both as 3-D tables and as projected spectra. A corresponding filter, which is 
necessary for comparisons with theoretical results and different from filters used for other 
data sets, is also available. 

Early in the second generation of the DLS program, measurements showed a previously 
unrecognized trigger inefficiency due to instantaneous rates much larger than the well con- 
trolled average rates. The microscopic duty factor of the beam, the instantaneous rates, 
and inefficiency varied on a time scale of hours. The efficiency could be reliably monitored 
only by the yield of the pairs themselves. A re-analysis of the first generation 4.9 GeV 
p+Be data found a correction factor of ~5 for the rate-dependent losses JRJ. Because of 
the unstable nature of this problem, it is unwarranted to assume that other first generation 
DLS data have the same correction factor. Although they do show evidence for significant 
rate-dependent inefficiencies, they lack sufficient information for calculation of appropriate 
correction factors. We suggest that the first generation data no longer be used for com- 
parison with theory. In the second generation data, the rate-dependent losses were greatly 
reduced by improvements to the electronics and the beam monitoring systems. For each 
data set, we measured the rate dependence of the pair yield to permit extrapolation to zero 
rate. After these improvements, a measurement of the differential cross section for p+p 
elastic scattering at 1.27 GeV (made concurrently with acquisition of p+p dielectron data) 
produced results in agreement with the known cross section []17|], and gave rate-dependent 
inefficiencies of < 15%. In the dielectron data reported here, run-by-run corrections for such 
losses ranged between 10-45%. 

Table I lists the data sets with the corresponding beam energies, pair statistics and abso- 



lute normalization uncertainties. The uncertainties are dominated by variations that result 
from using several methods to calculate the corrections for the rate-dependent efficiency. 
Cross sections were calculated from a single correction method, thus when comparing two 
data sets, the appropriate relative uncertainty is ~10-15%. 

The mass dependence of the differential cross sections is shown in Fig. 1. The systematic 
errors that are relevant to the shape of the spectra are displayed in this figure as added 
linearly to the statistical errors. These point-by-point systematics are independent of the 
normalization errors given in Table I. They are obtained from studies of the combinatoric 
background and studies of those acceptance corrections that were made in regions of phase 
space where the acceptance changes rapidly. Other effects, such as hadron contamination 
of the e + e~ sample, are negligible compared to the systematics shown. For masses below 
0.2 GeV/c 2 , an acceptance calculation with small bin widths has changed the shape of the 
data. (Henceforward we refer always to results within the DLS acceptance.) 

Fig. la displays the present data for the Ca+Ca cross section and the results of several 
calculations. The dotted line is from a BUU model || and shows the general trend of models 
that adequately represent our previous Ca+Ca data. For the present data the integrated 
cross section (M > 0.2 GeV/c 2 ) is ~7 times larger than both the earlier data and the model 
results. We attribute the difference between our two measurements to the uncorrected 
trigger inefficiency of the first generation DLS data. The agreement between models and 
the present data remains reasonable above 0.6 GeV/c 2 . In the mass range 0.2-0.4 GeV/c 2 , 
however, the models predict pair yields that are dominated by A and rj Dalitz decays but 
are significantly lower than our measurement. The shape of the yield is nonetheless quite 
similar to that of the 77-component of the BUU model shown in the dashed curve. 

To emphasize the information contained in the spectrum shape, we have used a simple 
model of ir° and r] production to calculate the mass spectra of the Dalitz decays. (Details of 
the model are discussed below.) The shapes of these spectra are insensitive to the parameters 
of the model, and the shape of the rj spectrum is in good agreement with that of the BUU 
model. Fitting the Ca+Ca data (0.05 GeV/c 2 < M < 0.375 GeV/c 2 ) with adjustable 
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amplitudes of these two shapes yields the dashed curve shown in the figure, with P(xD— 5%. 
The low probability of this fit may be due to the contributions from other sources. The 
overall agreement between curve and data shows that in this mass range the summed yield 
of pairs from those sources is either similar in shape to that from the ir° and r\ mesons, or 
relatively small, or slowly varying. 

As shown by the solid curves in the other panels of Fig. 1, the other data sets also are 
well represented (for 0.05 GeV/c 2 < M < ~0.4 GeV/c 2 ) by the same fitting procedure. 
We find that the fitted amplitude for each component, as well as the ratio of the cross 
sections, scale as the product of the projectile and target nucleon numbers, Ap-Ax- For 
the four reactions, the integrated cross sections (M < 0.35 GeV/c 2 ) scale as (Ap-Ay) Q with 
a = 1.06±0.01±0.02 (sys). 

A more direct comparison between the Ca+Ca and C+C data was obtained using the 
ratio of the cross sections as a function of the pair mass. This ratio, shown in Fig. 2, 
reveals two striking features. The first is that the ratio is independent of pair-mass for 
M < 0.4GeV/c 2 . Fitting the ratio to da/dM oc (A P - A T ) a , gives a = 1.01 ± 0.03+ 
0.04 (sys) - indicated by the line in the figure. This value is not inconsistent with the 
calculations of Ref. |7| for pairs from r\ decay [a = 0.87 ± 0.1) and A decay (a = 0.95 ± 0.1) 
produced in symmetric reactions ranging from Ca+Ca to Au+Au, and similar behavior has 
also been found for sub-threshold K + production at 1.0 A-GeV ||23| . Such a large value of a 
is unexpected for pairs from 7r° Dalitz decays, which are expected to dominate the low mass 
region M < 0.2GeV/c 2 . Because of the shape of the DLS acceptance, however, these low 
mass pairs are concentrated at rapidities > z/beam, and this may produce their strong Ap-Ap 
dependence. 

The second feature is the increase in the ratio for M > 0.5 GeV/c 2 , where fitting gives 
a = 1.40±0.13± 0.04 (sys). The data in Fig. 1 suggest that this value of a is due to high-mass 
contributions to the Ca+Ca data in addition to those producing the r^-like shapes seen in the 
other systems. Although some models have suggested that these pairs are primarily from 
7T + 7r - annihilations, Ref. |7j finds that the annihilation process should have a = 1.00+0.05 
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significantly different from our observation. Calculations of density-dependent effects of the 
medium on the pion propagator [I(J show a structure at ~0.5 GeV/c 2 in the p mass spectrum 



which could be contributing these additional pairs. Note that although these additional 
contributions could be from a source that produces pairs only with M> 0.5 GeV/c 2 , the 
present data do not rule out a source which contributes over a larger mass range but is 
visible only where other contributions (e.g. rj Dalitz decay) are small. 

Direct comparisons of mass spectra for other combinations of data sets are not made 
because at higher pair masses the rapidity distributions of the asymmetric systems are 
different from each other and from that of the symmetric systems. In the low-mass range 
these differences are not observed. 

Other data relevant to rj production in nucleus-nucleus collisions comes from measure- 
ments by the TAPS Collaboration of inclusive differential cross sections for 7r° and 77 pro- 
duction in 40 Ar+ nat Ca (and other) reactions at 1.0 A-GeV [fL8,19[|. The data span only a 



small rapidity interval near y cm . Meson production data in this limited rapidity range alone 
can not be used to calculate the Dalitz-decay yields of these mesons in the DLS acceptance 
because of the incomplete kinematic information. Thus there can be no model-independent 
statement regarding the fraction of the DLS pair yield that each of these mesons produces. 
The TAPS group calculated their total cross sections from a model where mesons were emit- 
ted isotropically from a thermal source at y cm . They obtained the temperatures of these 
sources by fitting their measured M r spectra. Both the magnitude and M t dependence of 
these data are reproduced by the BUU model ||. 

We have used the TAPS mid-rapidity data, inferred temperatures, and thermal model 
to calculate the contributions of their measurements within the DLS acceptance. This is 
the simple model referred to above. (We use a value of the 77 cross section 30% larger than 



the TAPS value because our data was taken at a higher beam energy [0.) The results 
are the histograms of Fig. la. The contribution from 7r°-Dalitz accounts for most of the 
Oct - \~ CcL cross section below 0.15 GeV/c 2 . The 77-Dalitz result, which is consistent with the 
77 component of the BUU calculation [g], is approximately 10% of the DLS yield near 0.25 
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GeV/c 2 . 

Significant (1 + a cos 2 ^) anisotropy is observed in both charged pion data pi|P2| , and 
sub-threshold K + production |23| . Because the DLS acceptance is small at y cm and peaks at 
2/beam, this angular distribution for meson production gives larger Dalitz decays yields than 
one which is isotropic. Our calculations show that the rj yield increases more than that from 
the 7T°. Without direct information on the angular distribution of rj production we cannot 
quantify the size of this effect. Nevertheless, given the significant uncertainties in both data 
sets ||24j| , we find that the TAPS data are not inconsistent with r] Dalitz decay contributing 
as much as 50% of the DLS pair production in the mass range 0.2-0.4 GeV/c 2 . 

In conclusion, we have presented dielectron measurements from Ca+Ca, C+C, He+Ca, 
and d+Ca collisions at 1.0 A-GeV. We find that the Ca+Ca cross section system is sig- 
nificantly larger than the original DLS measurement and the model calculations based on 
it. The low-mass cross sections from the four data sets reveal a mass- independent scaling 
of da/dM oc A p - A t suggesting similar dynamics for the dominant source mechanisms in 
reactions ranging from d+Ca to Ca+Ca. The shape of the low-mass spectra can be ap- 
proximated with the estimated pair distributions from tt° and r\ Dalitz decays. At higher 
pair mass, the ratio of Ca+Ca to C+C cross sections is much larger than the A p - A t ratio, 
indicating that a density-dependent mechanism(s) may be exhibited in this mass region. 
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TABLES 

TABLE I. DLS dielectron data sets. "Comb" denotes the measured combinatoric background, 
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40 Ca+ nat Ca 


1.04 


12800 


8102 


4698 ± 145 


± 30% 


12£j_|_natQ 


1.04 


4760 


1919 


2841 ± 82 


±30% 


4 He+ nat Ca 


1.04 


1929 


487 


1442 ± 49 


±30% 


2 H+ nat Ca 


1.06 


1828 


308 


1520 ± 43 


±40% 
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FIGURES 

FIG. 1. The DLS measurements of the dielectron cross sections from (a) Ca+Ca, (b) C+C, 
(c) He+Ca, and (d) d+Ca reactions. Panel (a) also contains the calculated signal from the BUU 
simulations of Ref. [9] (dotted line) and histograms showing the 7r° and rj decay contributions 
as estimated from the TAPS measurements and an isotropic thermal model. The dashed line 
represents the r\ component of the BUU calculation. The solid lines in all four panels show our fit to 
the low-mass data using the 7r° and r\ decay estimates (histograms) with adjustable normalizations. 

FIG. 2. The ratio of the cross sections from Ca+Ca and C+C collisions as a function of pair 
mass. The line on the plot is the fitted value of a from an assumed form da/dM oc(A p - At) a . The 
arrow indicates the kinematic limit for pair production in NN collisions. Asymmetric errors in the 
ratio occur where there are large fractional errors in the denominator; they are evaluated by Monte 
Carlo sampling. 
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Figure 1 . 
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